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We report the observation of the doubly Cabibbo-suppressed decay K^-k~h'^ in data from 

Fermilab charm hadroproduction experiment E791. With a signal of 59 ± 13 events we measured the 
ratio of the branching fraction for this mode to that of the Cabibbo-favored decay D®’ ^ A“7r®'7r'® 
to be B(D^ —> A®’7r“7r®') / B{D^ A“7r®'7r+) = (7.7 ± 1.7 ± 0.8) x 10“®. A Dalitz plot analysis 
was performed to search for resonant structures. 

The origin of the differences between the charm meson 
lifetimes is associated with their hadronic decays. While 
the semileptonic decay rates of the D® and are the 
same, the Cabibbo-favored (CF) hadronic decay rate of 
the D® is 3.2 times that of the D®". There are at least two 
possible sources for this difference. The CF D®" hadronic 
decay rate could be suppressed by destructive interfer¬ 
ence between spectator amplitudes containing indistin¬ 
guishable final state quarks. It is also possible that the 
CF D® decay rate is enhanced by non-spectator ampli¬ 
tudes which do not exist for the D®". 

For both hadronic CF D® decays and doubly Cabibbo- 
suppressed (DCS) D®" decays, all the final state quarks 
have different flavors, thus removing the possibility of 
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destructive interference. In the simplest picture, non¬ 
spectator amplitudes are small enough to ignore, and 
one would expect TDCsiD'^)/^CF{D^) ~ tan^^c and 
^dcs{D~^)/^Cf{D^) « 3.2 X tan'*6>c- These relations 
need not be satisfied if non-spectator amplitudes are also 
important. Doubly Cabibbo-suppressed decays can thus 
provide important insights into the D meson lifetime pat¬ 
tern. 

In this paper we report a measurement from Fermilab 
experiment E79I of the branching fraction for the DCS 
decay D+ ^ iF+7r“7r+. Throughout this paper, refer¬ 
ence to D+ and Df and to their decay modes imply also 
the corresponding charge-conjugate states. 

The data were recorded from 500 GeV/c tt~ inter¬ 
actions in five thin foils (one platinum, four diamond) 
separated by gaps of 1.34 to 1.39 cm. The experiment 
recorded 2 x 10^° events with a loose transverse energy 
trigger. 

The E791 spectrometer was an upgraded version of 
the apparatus used in Fermilab experiments E516, E691, 
and E769 Q. Position information for track and vertex 
reconstruction was provided by 23 silicon microstrip de¬ 
tectors (6 upstream of the target foils, 17 downstream), 
10 proportional wire chamber planes (8 upstream and 2 
downstream of the target) and 35 drift chamber planes. 
Momentum analysis was provided by two dipole mag¬ 
nets which bent particles in the horizontal plane. Particle 
identification was performed by two segmented threshold 
Cerenkov counters |^, allowing unambiguous identifica¬ 
tion of pions and kaons in the momentum range from 6 
to 40 GeV/c. 

After reconstruction, events with evidence of well- 
separated production (primary) and decay (secondary) 
vertices were retained for further analysis. The position 
resolutions along and transverse to the beam direction 
for the primary vertex were 350 ^m and 6 /rm, respec¬ 
tively. For 3-prong secondary vertices from decays, 
the transverse resolution was about 9 /im, nearly inde¬ 
pendent of the momentum; the longitudinal resolu¬ 
tion was about 360 ^m for a D'^ momentum of 70 GeV/c 
and increased roughly linearly with a slope of 30 ^m per 
10 GeV/c. 

We selected a generic Ktttt sample containing both 
DGS D+ ^ A'+7r“7r+ and CF D+ ^ A'“7r+7r+ decay 
candidates. The abundant CF decay was used to deter¬ 
mine the track and vertex selection criteria used in the 
search for the DCS decay. The criteria were chosen to 
maximize Ns/^/Nb, where Ns and Nb are the numbers 
of signal and background events in the sample. 

We required the secondary vertex to be well-separated 
from the primary vertex and located well outside the tar¬ 
get foils and other solid material, the momentum vector 
of the candidate D'^ to point back to the primary ver¬ 
tex, and the decay track candidates to pass closer to the 
secondary vertex than to the primary vertex. We used 
longitudinal separations normalized by their resolutions 
to reduce momentum-dependent effects. Specifically, a 3- 
prong secondary vertex had to be separated by at least 20 
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(Tl from the primary vertex and by at least 5 ctl from the 
closest material in the target foils, where the ctl are res¬ 
olutions in the measured longitudinal separations. The 
sum of the momentum vectors of the three tracks from 
this secondary vertex could not miss the primary vertex 
by more than 40/im in the plane perpendicular to the 
beam. We formed the ratio of each track’s smallest dis¬ 
tance from the secondary vertex to its smallest distance 
from the primary vertex, and required the product of 
these ratios for the three tracks to be less than 0.001. 

In addition to the selection criteria described above, 
we required Cerenkov particle identification for all three 
decay candidate tracks. The Cerenkov efficiencies and 
corresponding misidentification rates were measured us¬ 
ing the CF signal, in which the particle identification of 
the decay tracks can be determined from their charge. In 
this analysis the efficiency for correctly identifying kaons 
was 45%; the corresponding probability of misidentifiy- 
ing real pions as kaons was 2%. Since misidentification 
of the odd-charged pion candidate was a large source 
of contamination, we used a more stringent identifica¬ 
tion criterion for this track than that for the like-charged 
pion. For the odd-charged pion, the efficiency for correct 
identification was 57%, and the probability of misidenti- 
fying kaons as pions was 13%. For the like-charged pion 
the efficiency for correct identification was 85%; the cor¬ 
responding probability of misidentifying kaons as pions 
was 37%. The overall particle identification efficiency for 
£)+ ^ Ktytt was 22%. 

Due to particle misidentification and reconstruction er¬ 
rors, several other charm decays contributed to the back¬ 
ground in the K~*'7r~7r~^ sample. The major sources of 
charm background are listed below. 

a) 11+ and Dg decays with missing neutrals, such as 

Z1+ ^ Z1+ ^ iF“7r+7r+7r° and Dg (f)l'^u. 

In these cases 4-body decays produced 3-prong vertices. 
Monte Carlo simulations showed that, because of parti¬ 
cle misidentification, these events are spread smoothly 
across the entire it'+7r“7r+ mass spectrum. Their contri¬ 
bution was included with those of type (b) in a smooth 
background whose level was determined from the fit dis¬ 
cussed below. 

b) decays such as Ktt and ^ Kirmr. 

Such events passed the selection criteria for Ktttt when 
the reconstruction algorithm found two tracks from such 
a charm decay and combined them with a third track 
to form a spurious 3-prong vertex, or when one track 
was lost from a 4-body decay. False vertices created 
from it'“7r+ plus another track populate the 

iF+7r“7r+ mass spectrum above 2 GeV/c^, and were elim¬ 
inated by explicitly removing events with 77+71+ invari¬ 
ant mass between 1.828 GeV/c^ and 1.900 GeV/c^. A 
large fraction of the charm background originated from 

Ar“7r“7r+7r+ decays where one of the like-charged 
pions was lost and the remaining tracks were correctly 
identified. Since these events were concentrated below 
1.76 GeV/c^ the fit was restricted to the region above 
1.76 GeV/c^. Monte Garlo simulations indicated that 
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the remaining events from the —> K~Tr~'!T'^Tr'^ back¬ 
ground were smoothly spread across the mass 

spectrum. 

c) D'^ and Ds 3-body hadronic decays. These were 
the most problematic backgrounds because they produce 
structures in the it'+7r“7r+ mass distribution. Here, the 
3-prong candidates came from real 3-body decays whose 
reflections in the mass spectrum were con¬ 

centrated at shifted masses, except for the CF ^ 
K~Tr~^Tr~^ decay whose reflection was smoothly spread 
across the mass distribution. The Dg and 

£)+ ^ KKtt final states had very clean (j>TT components. 
We therefore eliminated events with K^K~ invariant 
mass in the range from 1.005 GeV/c^ to 1.035 GeV/c^. 
This removed 2% of true CF and DCS decays. 

The range of masses over which we could 

reliably model the charm background was 1.76 to 2.06 
GeV/c^. Within this interval, backgrounds of type a) 
and b) did not produce peaks. The structures resulting 
from background c) are shown in Figure 1. The param¬ 
eters for the reflection shapes were determined by inten¬ 
tionally misidentifying tracks in the background chan¬ 
nels. This was done with real data for Z1+ ^ it'“7r+7r+ 
and D^,Ds KKtt, and with Monte Carlo events for 
D+ TTTTTT. The net effect of these reflections is to pro¬ 
duce a small enhancement in the vicinity of the 

£)+ mass. 



FIG. 1. Distributions of the hadronic 3-body charm decay 
backgrounds that produce structures in the K'^'k~tt^ mass 
plot. The area under each curve represents the estimated 
number of background events of each type. 

The area under each curve in Figure 1 represents the 
estimated number of background events of each type, as 
fixed in the final fit. The amount of each background 
was determined by a rapidly-converging iterative proce¬ 
dure. The candidate events from the it'+7r“7r"'' sample 
were successively plotted as though they were K'^K~tt^ , 
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■K~^TT'^TT~ or K~, with an increasingly accurate de¬ 
scription of the feedthrough from the other channels. The 
result of this procedure was that in the K'^'k~'k^ mass 
range between 1.76 and 2.06 GeV/c^ there are 125 ± 
13 D'^ it'“7r+7r+ events, 25 ± 10 —> tt+tt+tt” 

events, 24 ± 11 Ug ^ K^K~'k^ events, and 15 ± 4 
Z1+ —!■ K^K~'k'^ events. 

The iir+7r“7r+ mass distribution for the final sample 
of decay candidates is shown in Figure 2. The spec¬ 
trum was fit to the sum of the reflections described 
above, a smooth function which describes the sum of 
all other backgrounds, a Gaussian function representing 
the D'^ 7f+7r“7r+ signal, and a Gaussian represent¬ 

ing the singly Gabibbo-suppressed Dg 7f+7r“7r+ sig¬ 
nal. The smooth background was modeled by an expo¬ 
nential function whose parameters were allowed to vary 
freely. The centroid and width of the Gaussian describ¬ 
ing the DGS signal were fixed to the values measured 
for the corresponding GF signal, 1.870 GeV/c^ and 12 
MeV/c^, as described below. The centroid and width of 
the Dg iF+7r“7r+ signal were fixed at 1.970 GeV/c^ 
and 12 MeV/c^. The number of D~^ iF+7r“7r+ signal 

events determined by the fit was 59 ± 13. 
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FIG. 2. Mass spectrum of candidate events with 

selection criteria optimized for the branching fraction mea¬ 
surement. The total background, including the contributions 
shown in Figure 1, is represented by the solid line outside 
the peaks and the dashed line under the peaks. The dotted 
histogram shows the mass distribution with tighter selection 
criteria. The events in the shaded area were used for the 
Dalitz plot analysis. 

The iF“7r+7r+ mass distribution for the sample of CF 
11+ ^ iF“7r+7r+ events, selected with the same criteria 
as the DGS signal and used as the normalization signal, 
is shown in Figure 3. The spectrum was fit to the sum 
of a linear background and a Gaussian signal whose pa¬ 
rameters were allowed to vary, yielding the values quoted 
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above. The number of K~'k'^'k^ signal events was 

found to be 7688 ± 90. 

Using Monte Carlo simulations with resonant and non¬ 
resonant decay modes, we found that the product of ac¬ 
ceptance and efficiency was the same for the CF and DCS 
samples within the ±2% statistical errors in the simu¬ 
lations. The ratio of branching fractions for the DCS 
D+ ^ and CF D+ ^ iF“7r+7r+ decay modes 

is, thus, given by the ratio of the measured numbers of 
DCS and CF signal events. 


B{D+ iF+7r“7r+) 
B{D+ K-TT+-K+) 


(7.7± 1.7±0.8) X 10“^ (1) 


The first error reported is statistical. The second is 
the systematic error, which was dominated by uncertain¬ 
ties in the background shape (8.5%), in the estimated 
number of background events used in the fit (4%), and 
the systematic error associated with particle identifica¬ 
tion (3.8%). The total fractional systematic error is 10%. 

Using the PDG value for the CF branching fraction, 
(9.1 ± 0.6)%, we find 

B{D+ iF+7r"7r+) = (7.0 ± 1.5 ± 0.9) x lO"'^, (2) 

where the fractional uncertainty in B{D^ 
has been added in quadrature with our systematic uncer¬ 
tainty in the ratio of branching fractions to determine the 
systematic error for the absolute DCS branching fraction. 

The value we have measured for the ratio of DCS to CF 
branching fractions is (3.0 ± 0.8)x tan'^dc, which agrees 
well with the simple spectator picture discussed in the 
introduction. It also agrees well with the recent result 
from Fermilab E687 0 for ratio (1), (7.2 ± 2.3 ± 1.7) x 
10 - 3 . 



FIG. 3. Mass spectrum of candidate D'^ —> 
events with final optimized selection criteria. The curve is 
a fit to a Gaussian signal and a linear background. The DCS 
D'*' ^ K^Tr~'K^ signal was normalized to this signal. 
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To study the DCS amplitudes which lead to the de¬ 
cay D'^ K^Tr~Tr'^, we have analyzed the Dalitz plot 
of a smaller but cleaner sample of events. Because so 
much of the background in the larger sample comes from 
misidentified charm decay, we used even more stringent 
particle identification on the odd-charged pion, which re¬ 
duced the particle identification efficiency from 22% to 
17%. We also required that each candidate’s proper de¬ 
cay time be greater than two lifetimes to suppress 
background from and Dg decays. We also explicitly 
removed candidates consistent with the mass hypothe¬ 
sis for D+ ^ This tighter selection gives 

a branching ratio consistent with that from the larger 
sample (equation (1)), for which selection criteria were 
chosen to maximize the projected sensitivity. The use of 
this cleaner sample, which contains 42 ± 9 signal events, 
reduced the systematic uncertainties from parametrizing 
backgrounds in the amplitude analysis. The Dalitz plot 
analysis was restricted to events found within 30 MeV jc? 
of the D+ mass (67 events altogether). These events cor¬ 
respond to the shaded area in the histogram at the bot¬ 
tom of Figure 2. The Dalitz plot of these events is shown 
in Figure 4. 



FIG. 4. Dalitz plot (M'^(JG^7r“) vs. M'^(7r^7r“)) of events 
found within 30 MeV/c'^ of the mass. This sample cor¬ 
responds to the events in the shaded region of the dotted 
histogram shown in Figure 2. 

The Dalitz plot was fit using the unbinned maximum- 
likelihood method. The D+ ^ TF+tt” 'Ti'^ decay ampli¬ 
tude was represented by a uniform nonresonant compo¬ 
nent plus relative amplitudes corresponding to the decays 
D+ ^ 7s:+p°(770) and £>+ ^ 7s:*°(890)7r+, 

2 

M = Anr + (3) 

i=i 

Each resonant component Rj was parameterized by 
a relativistic Breit-Wigner function, of constant width. 
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multiplied by a function describing the angular distribu¬ 
tion of decay particles. The nonresonant mode Ajvfl was 
chosen as the reference channel, fixing the scale for the 
relative fractions and the phase convention {Anh = 1). 
The fit parameters were, therefore, the relative phases 
6j and the real positive coefficients Qj for each resonant 
amplitude. The signal likelihood was obtained by multi¬ 
plying the ideal Dalitz plot density by a function describ¬ 
ing the geometrical acceptance and the reconstruction 
and event-selection efficiencies, including the removal of 
<j) K'^K~ and £)+ —> events described 

above. 

The background was represented by the sum of a con¬ 
stant term plus a product of two Gaussians in 
and 7r+7r“ masses accounting for the Ds K'^K~Tr~^ 
reflection, which is concentrated in the upper part of the 
Dalitz plot. The shape of the Dg reflection was obtained 
from Monte Carlo Dg events that passed 

through the same selection criteria as for the it'+7r“7r+ 
sample. We estimate 5 zt 3 Dg decays in the sample 
shown in Figure 4. 

The fit results are shown in Table 1. The decay frac¬ 
tions were computed by integrating the squared ampli¬ 
tude of each mode over the phase space, and then di¬ 
viding it by the integral of the square of the sum of all 
amplitudes. As shown in Table 1, the contributions of 
the three components are comparable. The two resonant 
modes are approximately in phase with each other and 
roughly 90° out of phase with the nonresonant part. 

Because of the Dg contribution to the cluster of events 
at the top of the Dalitz plot, the fractions corresponding 
to the A'“*’p°(770) and nonresonant modes are anticor¬ 
related with the estimated number of Dg K~^K~Tr~^ 
background events. Changing the estimated Dg back¬ 
ground to 8 events causes the K'^p°{770) fraction to de¬ 
crease by about 0.5 a, where a is the statistical uncer¬ 
tainty from the fit, and the nonresonant fraction to in¬ 
crease by the same amount. If, instead, the expected Dg 
background were 2 events, the K'^p'^{770) fraction would 
increase and nonresonant fraction decrease by 0.5 a. 
These uncertainties dominate the systematic error. 

We also attempted to include a D+ ^ Ar*(1430)7r+ 
amplitude in the fit, but this broad (T = 287 MeV) spin 
zero resonance could not be distinguished from the non- 
resonant amplitude. 

TABLE I. Results from the D'^ Dalitz plot 

fit. There is an anticorrelated systematic error of ± 0.07 on 
the fractions of the nonresonant and (770) modes, due 

to uncertainty in the estimated number of Dg —> K'^K~'k^ 
background events. This uncertainty in Dg background has 
negligible effect on the phases. 


mode 

phase (radians) 

fraction 

nonresonant 

0 (fixed) 

0.36 ±0.14 ±0.07 

A*°(890)7r+ 

1.8 ±0.5 

0.35 ±0.14 ±0.01 

K+p° (770) 

2.0 ±0.4 

0.37 ±0.14 ±0.07 



(GeVVc‘) 


The fit results are also shown in Figure 5 as Monte 
Carlo density plots. On the left is the generated signal 
probability distribution function (pdf); on the right is 
the signal plus background pdf, now multiplied by effi¬ 
ciency and acceptance functions. Due to the more strin¬ 
gent particle-identification criterion adopted for the odd- 
charged pion candidate, the overall efficiency was reduced 
at low values of the squared masses. This can be seen 
clearly by comparing plots a) and b) in Figure 5. 

A Monte Carlo technique was used both to test the fit¬ 
ting procedure and to assess the goodness-of-fit. A large 
number of Ar''"7r“7r“'' samples was generated according to 
the overall probability distribution function, using as in¬ 
put parameters the phases and coefficients given by the 
fit to the real data. Each of these Monte Carlo samples 
was then fitted, and the distributions of the resulting fit 
parameters were plotted. The average values of all fit 
parameters were the same as their input values, and the 
rms spread in each fit parameter distribution was close to 
the error on this parameter given by the fit of real data. 
The fraction of Monte Carlo samples for which the value 
of —2 In {Lmax) exceeds that of real data, where Cmax is 
the maximum value of the sample likelihood, estimates 
the confidence level of our fit to be 52%. 


(a) (b) 



FIG. 5. Monte Carlo simulations of the 
Dalitz plot for (a) generated signal and (b) accepted signal 
plus background. Both use the signal parameters of the fit 
reported in Table I. Plot (b) also includes the background 
measured for the data of Figure 4. In these plots, the lin¬ 
ear dimension of the boxes is proportional to the number of 
events. 

In spite of the large errors on the fractions, the pattern 
shown in the DCS D+ ^ Ar+7r“7r+ decay appears to 
differ qualitatively from the corresponding CF decay, in 
which the nonresonant component corresponds to about 
95% of the branching fraction . 

In summary, Fermilab experiment E791 has measured 
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the ratio of branching fractions for the doubly Cabibbo- 
suppressed Z?'*" ^ and the Cabibbo-favored 

ZZ+ ^ Zf“7r+7r+ decay modes, B{D'^ Zf+7r“7r+) / 

B{D+ K-7T+7r+) = (7.7 ± 1.7 ± 0.8) x IQ-^, corre¬ 
sponding to (3.0±0.8) tan^ 0c- A Dalitz plot analysis in¬ 
dicates that the DCS signal is composed of approximately 
equal amounts of Z1+ ^ it'+p°(770), 11+ ^ Zf*°(890)7r+ 
and nonresonant modes. Using the measured fractions 
from Table 1 and the 11+ ^ Zf+7r“7r+ branching fraction 
from equation (2), we obtain B{D^ itr*°(890)7r+) = 
(2.5 ± 1.2) X 10-*, B{D+ K+p°{770j) = (2.6 ± 
1.3) X 10“*, and B{D^ —> nonresonant Zf+7r“7r+) = 
(2.5 ± 1.3) X 10“*; statistical and systematic errors have 
been added in quadrature. 
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